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ABSTRACT 

At radio wavelengths, scattering in the interstellar medium distorts the appearance of astronomical 
sources. Averaged over a scattering ensemble, the result is a blurred image of the source. However, 
Narayan & Goodman (1989) and Goodman & Narayan (1989) showed that for an incomplete average, 
scattering introduces refractive substructure in the image of a point source that is both persistent 
and wideband. We show that this substructure is quenched but not smoothed by an extended source. 

As a result, when the scatter-broadening is comparable to or exceeds the unscattered source size, the 
scattering can introduce spurious compact features into images. In addition, we derive efficient strate¬ 
gies to numerically compute realistic scattered images, and we present characteristic examples from 
simulations. Our results show that refractive substructure is an important consideration for ongoing 
missions at the highest angular resolutions, and we discuss specific implications for RadioAstron and 
the Event Horizon Telescope. 

Keywords: radio continuum: ISM - scattering - ISM: structure - Galaxy: nucleus - techniques: 
interferometric 


1. INTRODUCTION 

Radio-wave scattering in the turbulent interstellar medium (ISM) produces familiar effects: scintillation in frequency, 
time, and position. This scintillation has two distinct branches: diffractive and refractive. Diffractive scintillation is 
narrowband, short-lived, and is quenched by a source exceeding the diffractive scale, corresponding to the resolution 
of the scattering disk when viewed as a lens. As a result, diffractive scintillation is only observed for pulsars, masers, 
and a few extremely compact quasars. In contrast, refractive scintillation is wideband, persistent, bu t is quenched 
only w hen the angular size of the source exceeds that of the scattering disk. Reviews by fRickett (1990) and Narayan 
( 1992) outline scattering theor y and observations in both of these regimes. 


Narayan & Goodman (1989 hereafter NG89) and Goodman & Narayan (1989 hereafter GN89) uncovered a sur- 
prismg refractive effect: substructure in the scattered image of a point source. Specifically, they showed that this 
substructure contributes noise to interferometric visibilities on baselines long enough to resolve the smooth, ensemble- 
average scattered image. This noise is wideband and p ersists over the refra ctive timescale (i.e., the time for the 
scattering material to move across the scattered image). Gwinn et al. (2014) recently provided a dramatic demon¬ 
stration of this effect through the discovery of substructure in the heavify scattered image of the Galactic center 
supermassive black hole, Sgr A*, at 1.3-cm wavelength. 

We now extend the theoretical framework that describes refractive substructure to accommodate short interferometric 
baselines, anisotropic scattering, and extended source structure. We show that, surprisingly, substructure in the 
scattered image can occur at much finer angular scales than those of the unscattered source. Thus, although an 
extended source smoothes the diffraction pattern in the observing plane according to the familiar convolution action of 
scattering, it does not smooth the scattered image but merely reduces the depth of fluctuations. As a result, refractive 
substructure can remain an important consideration for observations in which the scattering is somewhat subdominant 
to intrinsic structure, and it can introduce spurious compact features into resolved images of extended sources. 

We begin, in (j^by reviewing some basic principles of scattering and scintillation. Next, in ^ we explore interfero¬ 
metric visibilities in different averaging regimes. In ^ we define the scattered image and consider how its appearance 
is affected by properties of the source and the scattering. We also derive expressions that allow for efficient numerical 
computation of scattered images and provide characteristic results from simulations. In we derive specific observ¬ 
able properties of the refractive noise - flux modulation, image wander, and substructure in scattered images - using 
the second moment of the interferometric visibility modulus (derived in Appendix . In ® we consider implications 
of refractive noise for two specific missions: RadioAstron and the Event Horizon Telescope. We summarize our findings 
in ^ 

2. THEORETIGAL BACKGROUND 
2.1. Interstellar Scattering and Scintillation 

Scattering of radio waves in the interstellar plasma arises from small-scale fluctuations in electron density. The 
resulting variations in refractive index produce variations in phase of the electromagnetic wave. As a result, each 
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scattered ray has both geometrical and stochastic contributions to phase. 

In many cases, the scattering can be well-approximated by a single thin phase-changing screen. The change in 
phase by the screen is (^(x), where x is a transverse coordinate on the screen. The statistical characteristics of the 
scattering and scintillation can then be related to statistical characteristics of the phase fluctuations, either through 

a spatial structure function D^{x) = /[</) (xp -I- x) — (/)(xo)]^ ) or, equivalently, through the power spectrum Q(q) of 


phase fluctuations. _ 

A variety of evidence (e.g., Armstrong et al.|1995 ) suggests that these functions are well approximated as power laws 
- D^{yi) oc |x|“ and Q(q) oc |q|~^“~'~^'* - spanning the immense range of scales from ~10^ km to >1 Auj^On sho rter 
scales, the phase fluctuations are smooth, so D^{yi) oc |x|^, while on longer scales, is constant (TatarS;ii||l971 ). 
The density fluctuations, and their power-law form, m ay reflect the effects of a cascade of Alfven-wave turbulence 


(Goldreich & Sridhar 1995 Lithwick & Goldreich 


|2001 ). Indeed, scattering often follows the Kolmogorov scaling 
in power-law index 0 = 5/3. In this picture, the cascade is initiated 
by driving forces at large spatial scales, the “outer scale” Tout, and is terminated by dissipation at a minimum scale. 


expected tor this and other types of turbulence, wit 


the “inner scale” rin. Some evidence su 

jjgests that the inner scale for interstellar scattering material is 100 to 300 km 

(| Spangler & Gwinn 

1990 

Rickett et al. 

|2009), although only a few, highly-scattered lines of sight have been studied. 

i'he inner scale anc 

turbulence properties may differ for paths along and perpendicular to the large-scale magnetic 


field. We provide simple expression for D^{x) and Q{q) that accommodate anisotropy and inner and outer scales in 
Appendix and 

If the power-law index o < 2 then the spectra is said to be “shallow,” whereas for a > 2 the spectra is said to be 
“steep.” In this paper, we restrict our attention to shallow spectra. 


2.2. Scintillation and Averaging Regimes 

In the strong scattering regime, there are three important length scales on the scattering screen. The phase coherence 
lengthy rg, corresponds to the separation between two points for which the root-mean-square (rms) screen phase 
difference is 1 radian: D^lrg) = 1. The Fresnel scale, rp, defines the lateral scale at which the change in geometrical 
phase relative to that of the direct path is 1/2 radian and is given explicitly in Eq. below. Finally, the refractive 
scale, tr, determines the size of the scattered image of a point source: tr = r|/ro. The strong scattering regime is 
defined by the condition rg rp r^i. 

These length scales delineate the two scintillation regimes, diffractive and refractive, introduced in ^ Diffractive 
scintillation is dominated by fluctuations on the scale of rg whereas refractive scintillation is dominated by fluctuations 
on the scale of tr. Thus, diffractive scintillation decorrelates over a fractional bandwidth of '^rg/rpi while refractive 
scintillation decorrelates over a fractional bandwidth of unity. Likewise, if u_l denotes the characteristic transverse 
velocity of the scattering material, then rg/v± gives the diffractive timescale while rpi/v± gives the refractive timescale. 

Following the treatment and nomenclature of |GN89| we will consider three types of averages for quantities such as 
interferometric visibilities (denoted 14) and images (denoted /j,). A snapshot quantity (e.g., 14s) averages over source 
and background noise for a single realization of the scintillation pattern. An average quantity (e.g., 14) averages also 
over diffractive scintillation but not refractive scintillation. An ensemble-average quantity (e.g., 14a) averages over 
both diffractive and refractive scintillation. We will also use these subscripts to denote respective averages (e.g. (.. .)ss 
represents a snapshot average). 

2.3. The Scalar Electric Field 


To derive observable consequences of scattering, we must first determine the scalar electric field tp{h) at a transverse 
coordinate b in the observing plane. Using the Fresnel diffraction integral, this field can be written as (see, for instance, 
Gwinn et al.||1998 ) 

^(b) = -^/ d2xe*[(2fe)|b->^l'+<^W] / = 

./screen ^src 


In this expression, a large, constant phase is absorbed into the stochastic source field ^src(s), and we have chosen 
the amplitude to simplify later calculations. Here, (/'(x) is the screen phase, D is the characteristic Earth-scatterer 
distance, R is the characteristic source-scatterer distance, k = 2 tt/X is the wavenumber, and rp is the Fresnel scale, 
defined bjj^ 


rp = 


DR 1 
D Rk 


( 2 ) 


We will write most subsequent equations in terms of rp and the effective magnification M = D/R of the scattering 
screen. 


3. INTERFEROMETRIC VISIBILITIES IN DIFFERENT AVERAGING REGIMES 


We now consider the behavior of interferometric visibility (i.e., the electric field covariance) in the various averaging 


regimes outlined in (2.2 


^ Note that a is used inconsistently throughout the literature, with some authors choosing a as the index for Q(q) so that oc |x|“ ^ 

^ Note that some authors define the Fresnel scale differently, often with the substitution k — ^ fc/(27r). 
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Figure 1. Geometry relating the source, scattering, and observer. The scales rgrc rjmg correspond to the coherence lengths in the 
observing plane of the electric field from the source and ensemble-average scattered images, respectively. These scales are related to the 
angular sizes of these images via rsrc ~ O.37A/0src and rimg ~ O.37A/0img, with rsrc > ^mg- Also, ro denotes the transverse scale on the 
scattering screen over which the rms phase difference is one radian; the scattered angular size of a point source is 0scatt ~ 0.37A/ [(1 M)ro]. 


3.1. The Snapshot Visibility 

The simplest averaging regime is the snapshot - an average over noise of the source for a fixed realization of 
the scattering. We assume that the source is spatially incoherent with an intensity /src{s): (V;src(s)'i/'src(®0)ss = 
/src(s)(5 (s — s'). Then, the snapshot visibility on a vector baseline b centered on bg is (see Figure]^ 

y,(b;bo) = {^{ho - h/2)r{ho + b/ 2 ))^^ (3) 

= -f— [ + f d2g (x^-xi ).s 

dTT^r^ J J 

= e*5’'F'[(^?-^2) + THlW'(xi+X2)-2^.(b+Xi-X2)]gZ[,^(xi)-0(x.)]y^^^((j 

J 


In the final expression, we have applied the Van Cittert-Zernike Theorem (e.g., Thompson et al.||2001 
remaining source integral (a Fourier transform of the source brightness distribution) with its equiva. 
interferometric visibility I4rc corresponding to the unscattered source image: 


Kirc(b) oc J d2s4rc(s)e^’''(o+«) (^). 


) to replace the 
ent form as an 


(4) 


Because our later results are all expressed as fractional quantities, overall normalization is insignificant. To simplify 
calculations, we will adopt the normalization I4rc(0) = 1- 
Note that the intensity of the diffraction (or illumination) pattern in the observing plane takes the form 


4s(bo) = (i/’(bo)i/’*(bo)) 
_ 1 
Att'^ 


(5) 


^ [ d2x^d2^2e'5'’p'[(^i"^')”^T5^4xi-X2)]gd0(xi)-0(x2)] f rf2ggi^(x2-xi).sj^^ /gX 

J J 

= J X J d2xid2x2e'5’'F'[(^?-^2)]ed0(xi)-0(x2)]g*fc(* + T)(x2-xi)| _ 

The term in braces is simply the point source diffraction pattern at a transverse coordinate bp + Ms. So, for an 
extended source, we recover the well-known result that the diffraction patter n in the observing plane is the con volution 
of the point-source diffraction pattern with a magnified image of the source ( SalpeterjlQB? Cohen et al. 1967). Hence, 
an extended source will smooth out any features in the diffraction pattern that are hner than the scales of the inagnified 
source image. However, in Qwe will demonstrate that an extended source does not smooth the scattered image. 


3.2. The Average Visibility 

The snapshot visibility, corresponding to Eq.[^ represents the cumu lative e ffect of three contributions: the ensemble- 
average visibility, diffractive noise, and refractive noise. As shown by|GN89 and illustrated in Figure the diffractive 
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contribution to the integral comes from regions with |xi — X 2 I 3> Tq (see Appendix 1 ); the ensemble average and 
refractive noise arise from the region with |xi — X 2 I < tq. Indeed, this conclusion is apparent because an extended 
source limits the con tribution from widely-separated pairs of points and thereby quenches the diffractive scintillation 
(Gwinn et al. 1998). Thus, a extended source with size that is significantly larger than the diffractive scale will 
immediately take shapshot visibilities into the regime of the average visibilities. When the angular size of the source 
exceeds the refractive scale (i.e., when the scattering becomes subdominant to source structure), the source will impose 
an even tighter restriction |xi — X 2 I < (1 -|- M)~^r < tq that reflects the shorter correlation length of the source 
electric field at the screen. We define rimg more precisely in ^ 

A coarse examination of Eq. reveals most salient propmies of the noise in average visibilities. For instance, 
we can readily understand why average visibilities are quite insensitive to integration in frequency or time - it is 

For example, any significant effect from a shift 

> 1. Thus, short baselines must be 


a consequence of the restriction |xi — X 2 I (1 -I- M) ^nmg < ?’o- 
of the baseline center Abp requires (1 -|- |Abo| -I- (1 -I- M)~^r. 


- 2„2 


displaced by ^(1 -I- M)^rp/rii„g, which is the appropriately magnified size of the scattered image. Long baselines 
must be displaced by an amount that is smaller by ~ (1 -I- M)~^ri^g/b. Because baselines with b ro resolve the 
ensemble-average scattering disk and, therefore, see little flux, the former is the more typical circumstance. Thus, to 
influence accessible refractive metrics, baseline centers must be displaced by the refractive scale. An extended source 
increases the spatial correlation scale, reinforcing the notion of a smooth diffractive pattern in the observing plane, 
the convolution of the source intensity distribution with the point-source response, as discussed in §3.1[ 

We can likewise consider the effects of averaging in frequency. To incorporate frequency dependence, there are two 


necessary modifications: rp oc V~X and (j){x) oc A. The additional phase from a wavelength change AA is then given 
W (^1 ~ ^2) + + M)~^h ■ (xi -I- X 2 ) -I- [(/'(xi) — (()(x 2 )]}. In the refractive regime, |xi — X 2 I < ro, so 

(j){xi) — (()(x 2 ) < 1. Hence, the phase of the first term is the second is and the third is Thus, 

the fractional bandwidth of refractive noise is of order unity, except on long baselines where it becomes ^Vimg/b. Note 
that an extended source will decrease the fractional bandwidth of refractive noise on long baselines but will not change 
the fractional bandwidth of refractive noise on short baselines. 

However, while the baseline center bp must shift by the refractive scale to incur a significant change in the average 
visibility, the vector baseline b only needs to change by a distance of ^rimg. Thus, an extended source will decrease 
the coherence length of the refractive noise in the visibility domain. This property already suggests the presence of 
intense substructure in the scattered image, which we explore in ^ 

Because the current paper emphasizes refractive effects, we will simply set: 


bo = 0 


( 6 ) 


and ignore frequency dependence for the remainder. 


3.3. The Ensemble-Average Visibility 

To obtain the ensemble-average visibility, we must average over many realizations of the scattering screen, defined 
by the Gaussian random field (^(x). To calculate this average, the characteristic function of a zero-mean Gaussian 

random variable provides a convenient identity: = g -2 ^ 

then change the integration variables to y = Xi — X 2 and x = (xi -|- X 2)/2 to give 

Fea(b) = f (1 + M)y). (7) 

F t/ 

Integration over x is then trivial, using the identity J d^xe^^'^ = (27r)^(5^(y). Thus, because D^{y) = D^{—y), 

Fea(b) = 


( 8 ) 
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Hence, we recover the well-known result that the ensemble-average visibility for a scattered source is the product of 
the point source visibility and the unscattered source visibility. A familiar corollary is that, in the ensemble-average 
regime, scattering convolves the unscattered source image with a smooth scattering kernel. 


Indeed, because the ensemble-average scattering is deterministic and has a positive kernel, it is readily 

inverted when known (Fish et al. 2014). Hence, in the ensemble-average regime, one can reconstruct the unscattered 
image of the source by dividing nieasured visibilities by the scattering kernel. With this approach, the residual effects 
of scattering are simply to amplify the thermal noise on long baselines. 

Moreove r the conv olution action of scat tering leads to other convenient properties in this regime, as discussed by 
Fish et al. (2014) and Johnson et al. (|2014 ). For instance, because the kernel is real and positive, visibility phases are 


unaffected by scattering in the ensemble-average regime, even if the scattering kernel is not known. Also, quotients of 
visibilities on equal baselines will be unaffected, such as fractional polarization in the visibility domain (because the 
scattering is not significantly birefringent). 


4. THE SCATTERED IMAGE 
4.1. Definition of the Scattered Image 

We now consider the appearance of the scattered image. However, there are two immediate difficulties in precisely 
defining the scattered image. The first is that it may depend so sensitively on observing position that any aperture 
sufficient to resolve the image will necessarily span different image elements. The second is that the scattered image 
is not incoherent, and so is not trivially related to visibilities via the Van Cittert-Zernike theorem. 

The first of these concerns is problematic for snapshot visibilities. However, we have seen that average visibilities are 
insensitive to changes of observing position (i.e., baseline center bg) that are less than the refractive scale. Thus, the 
average image is well-defined in this regard. For the second concern, because the coherence length on the scattering 
screen is ~ rg tr, the scattered image can be considered effectively incoherent for all baselines of length |b| ^ 
XD/ro ~ tr. 

Hence, we can define the average image according to the Van Cittert-Zernike theorem applied to average visibilities: 

/avg(x)cx y'd2bK,g(b)e-2”(5)-(S) = y’d2bVavg(b)e”^"^ ". (9) 


Note that we have calculated the image at a distance of the scattering screen, D, rather than at a distance of the 


ING891 

;he snai 


and allows linear coordinates to be easily compared on the screen, 
snapshot visibility on bg, in accord with the discussion of H3.2| (a shift 


source. This convention was also adopted by 
Note also that we ignore the dependence of t 
in bg is equivalent to a shift of the scattered image). Although this image is rather different from that produced by a 
physically realizable interferometer or aperture, it can be filtered or averaged to reproduce their behavior flexibly and 
easily. 

Observe that the integrand of Eq.|^is conjugated when b — >• —b, so Javg(x) G M, although the image is not necessarily 
positive. 


4.2. The Scattered Image of an Extended Source 

We now derive an expression for the scattered image of an extended source in the average-image regime. Here, we 
require only that the source quenches the diffractive scintillation so that the snapshot- and average-image visibilities 
are equal (see The preceding definition for the scattered image is especially useful in this case because the snapshot 
visibility has a straightforward dependence on baseline. After substituting 14s (b) (given by Eq.[^ for 14vg(b) in Eq.[^ 
the integral over b gives a delta function with argument proportional to x — (xi -|- X2)/2. Thus, the only cont ribut ion 
to the snapshot image at a location x is from pairs of points on the screen that are centered on x. As in §3.3[ we 
can change to variables given by the average and difference of xi and X 2 . Integrating over the former leaves a single 
remaining integral over y = X 2 — xi: 


,(x) oc / Vsrc ((1 + M)y) e 






x+^y 


)], 




( 10 ) 


This form is especially convenient for numerical simulations. In §4.2.3^ we will derive an even simpler approximate 
representation, which eliminates the remaining integral. 

Note that if the screen-phase term exp {i (x — ly) — ^ (x -|- ly)] } depended only on y, then /avg(x) would be the 


Fourier transform of the product of a position-independent scattering kernel and the source visibility (compare Eq. 10 
with Eq. 1^ ; consequently, the scattered image would be the convolution of the point-source response with an image of 
the source. Perhaps unfortunately, the screen-phase term depends on x so the scattering does not act as a convolution 
for average images. However, after an ensemble-average over the screen phases, (exp {z [^(x — ^y) — (j){x + ^y)] = 

so that we again recover the expected convolution action of scattering in the ensemble-average regime. 


4.2.1. Existence and Persistence of Refractive Substructure 

The arguments of the preceding section raise the question: why should refractive substructure exist at all? If the 
screen phase, ffxf, decorrelates over the scale of rg at the screen, why does each small region of scale rg not produce an 
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Figure 3. Simulated images showing the effects of scattering for two sources: a circular Gaussian distribution {top) and a ring (bottom). 
The leftmost panel in each shows the unscattered source; the rightmost shows the ensemble-average scattered source. The three central 
panels show three different average images, each with a Kolmogorov index in the inertial range. The second and third have a large inner 
scale (rjn = 200ro ~ O.lGrp^) and a small outer scale (rout = 500ro ~ 0.41rR), respectively. Each image has rp = 35ro and M = 1 and was 
calculated using a scattering screen with 2^“^ x random phases, with the same random seed for both source mod els. W e used rin = 5ro 
in each case (except the large inner scale examples) so that the description and numerical reduction of described in |4.2.3|was appropriate. 
The color scale is linear and extends to the brightest pixel in each image; the angular size of each image is identical. Observe that the 
greater extent of the Gaussian than the ring r esults in a shallower ^nature of substructure in the average images but does not blur the 
substructure, in accord with the discussion of See Appendix [3 for additional details about the simulations. 


independent variation from the average image, so that any larger-scale structures are completely random in character? 
As we will now demonstrate, although the phase decoheres over the scale tq, the phase gradient remains correlated over 
much larger scales. This long-range correlation accounts for the existence and persistence of refractive substructure. 

Specifically, note that the autocovariance of the directional derivative of phase can be related to the second derivative 
of the phase autocovariance: 

([a:r(/>(xo)] (xo -b x)]) = -dl ((?!) (xo) (xo -b x)) = d‘^D^{x.). (11) 


Thus, for a power-law index a, the phase-slope autocovariance only falls as |x| (e.g., as 

spectrum). More p recisely, note that the phase gradient causes a shift of the image by ^r|V())(x) (we will show this 
explicitly in j 4.2.3| ). Thus, the autocovariance of the fractional image shift, relative to the refractive scale, of points 
offset by x is ~(|x| A large inner scale extends this correlation length, while an outer scale introduces a 

sharp cutoff for the correlation. The broad correlation of phase gradient leads to the coherence of substructure on 
scales much larger than rp and allows it to persist over relatively long temporal averages which, with the assumptions 
of a “frozen” phase screen, are equivalent to partial spatial averages. 


x| for a Kolmogorov 


4.2.2. The Role of Source Structure 

We now address another fundamental question: does an extended source smooth the scattered image? To answer 
this, note that the effect of an extended source in Eq.[^is to restrict the range of y to small values. For a given location 
X on the image, this restriction decreases the variations of exp {i [(() (x — ^y) — ^ (x -b ^y)] }, which are the source of 
the refractive substructure (and the scatter broadening). This phase difference quickly approaches 0 when the angular 
size of the source exceeds the ensemble-average angular size of a scattered point source, thereby quenching the refractive 
substructure. However, because the image coordinate x determines the center of the separation of contributing pairs 
of screen points in Eq. that produce the refractive substructure, the only role of source structure is to diminish 
the effects of this phasedifference at each x independently. Source structure will reduce the screen-induced phase 
fluctuations but will not spatially smooth them. 


4.2.3. Approximating the Scattered Image 

A limiting case provides considerable insight and a powerful tool for simplifying numerical work. Specihcally, 
suppose that the screen phase fluctuations are smooth over the range of co her ence for the source field at the screen: 


Ksrc ((1 + M)rin) ^ 1. Then, over the relevant range of integration in Eq. 10 (j){x + ky) — — hy) ~ y ' V(()(x). 
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Figure 4. Simulated images showing the effects of source size on a scattered image of a circular Gaussian source, as in Figure Four 
cases are shown: the leftmost has a source that has the same angular size as the ensemble-average scattered image of a point source, the 
next three are for identical scattering but an intrinsic source that is 2, 4, and 8 times larger. In each case, the scattering is Kolmogorov 
with a small inner scale and infinite outer scale, a Fresnel scale rp = 17.9ro, a magnification M = 1, and a scattering screen with 2^'^ x 2^“^ 
random phases. The color scale is linear and extends to the brightest pixel in each image; the angular size of each image is identical. 
Even when the scattering is highly subdominant to the intrinsic structure, the effects of scattering-induced substructure are still readily 
apparent. 


With this approximation, Eq. becomes 


4s(x) « J (fy Vsrc ((1 + M)y) 

= /src (x + r|V(/)(x)) . 


( 12 ) 


In this case, refraction “shuffles” the image, moving each image element to a place nearby: the scattered image 
brightness at a location x is given by the brightness of the un scattered source at x + rpV(()(x), in accordance with the 


brightness theorem of geometrical optics (Born k. Wolf 1980). Since |V(/)| ^ I/tq, the shuffling occurs over a region 
spanning the refractive scale, centered on x. Again, we emphasize that x is a transverse coordinate at the scattering 
screen (not the source) with a corresponding angular coordinate x/U. 

This approximation vividly illustrates the breakdown of the convolution action of scattering in the average image, 
as discussed above. Conversely, in the ensemble-average limit, the averaged “shuffling” behavior leads to a smooth 
position-independent blurring around each image coordinate to reproduce the familiar convolution of point-source 
response and image, now in the ensemble-average regime. 

Because the approximation of Eq. [^eliminates the remaining position-dependent integral, it enables rapid estimation 
of the effects of scattering on any image after generating an appropriate random phase screen. Figure shows some 
examples of scattered images that utilize this approximation. 

4.3. The Imprint of Scattering Characteristics on the Scattered Image 


The discussion of (4.2 allows us to asses s how the scattered image is modified by properties of the scattering (see 
Figure [^ for examples). For instance, in (4.2.1 we showed that the large-scale coherent structures in the scattered 
image are related to the long-range correlation of the phase gradient V^(x) across the image, which is proportional to 
the second derivative of 

Thus, if there is a finite outer scale, then the phase gradient must decohere over displacements of "^Cout to avoid 
excessive growth of the phase and divergence of the phase structure function. This decoherence quickly destroys the 
coherent features in the scattered image on scales larger than rout and eliminates power on baselines too short to resolve 
the scale of the decoherence (|b| ^ ~ f^rg). On the other hand, a large inner scale increases the correlation 

length of the phase gradient. This increase results in larger-scale coherent features in the scattered image, or increased 
“patchiness.” 

Similar considerations apply to the po wer -law index a. Because larger values of a < 2 have an increased correlation 
length in the phase gradient (from Eq. [Tl|), they produce a higher level of refractive noise on short baselines but a 
steeper fall in the noise with increasing baseline. 

5. MANIFESTATIONS OF REFRACTIVE NOISE 

We now discuss observable manifestations of the refractive noise, how each is quenched by an extended source, and 
how each reflects properties of the scattering. Our discussion relies on expressions for the variance of the snapshot 
visibility that are derived in Appendix [^ We outline how to obtain exact expressions numerically (using results in 
Appendix B.2| but primarily focus on developing intuitive understanding, approximations, and scaling relationships. 

In §5.1[ we discuss the most familiar refractive effect, modulation of t he to tal image flux. Next, in §5.2[ we derive 
the refractive contributions to image wander and distortion. Finally, in §5.3[ we derive expressions for the visib ility 


noise on long baselines, reflecting refractive substructure in the scattered image. We summarize these results in (5.4 
and in Table [T] 
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5.1. Refractive Flux Modulation 
Th e most familiar refractive manifestation, modulation of the total flux, has been studied in depth by many authors 


Shishov 

1974 Rickett et al. 

1984 

Goodman & Narayan 

1985 


__ __ 

the visibility measured by a zero- baseline mterterometer, we can use our formulation to estimate the flux modulation 

providing a valuable point-of-contact with previous work. 

To derive an explicit estimation the flux modulation, we will approximate the ensemble-average image by a Gaussian 
with a characteristic scale nmg in the visibility domain: 14a (b) = e^sd^l/orng) 


The reciprocal of rimg 

e„ ^ 


is then 
0.37^. 
and the 


proportional to the full width at half maximum (FWHM) of the scattered image, 6*img- ^^img 
In turn, ^img is approximately given by the quadrature sum of the FWHM of the unscattered source, 0src, 

FWHM of a point source, ^gcatt- 

With these substitutions, is it straightforward to evaluate Eq. on the zero-baseline. The rms flux variations are 
given by a/((A/)^) cx (6'scatt/6'img)^~“^^, in agreement with previous derivations (e.g.. 


N arayanj 1992 1) . 


From Eq. 


32 


and Eq. 34 it is evident that an inner scale only weakly affects the flux modulation until rjn > r^ro/rinig; 
the inner scale must be comparable to the transverse extent of the ensemble-average image to have a significant effect. 
However, the imprint of a small outer scale is much more severe - as th e outer scale becomes much less than tr, the 
rms flux modulation will be quenched as (see Eq. 35). 

5.2. Refractive Image Wander and Distortion 

Observations that are sensitive to the absolute image position will see refractive image wander, and observations on 
short baselines that begin to resolve the scattered image will see large-scale image distortion. In general, these e ffects 
will be correlated with each other and with the flux modulation (|Blandford & Narayan|1985 |Romani et al.|1986 ). We 
will now show that these effects are quit e gen erally related to retractive noise on short mterierometric baselines, and 
we will then use the results of Appendix |B. 2| to estimate their magnitude. 

Working in dimensionless baseline u = b/A for the visibility V (u) and angular coordinates for the source image 
I{0), the Van Cittert-Zernicke Theorem takes the form 


v{u) = j d^ei{e)( 


2'Kiu-9 


(13) 


For a short baseline - i.e., a baseline that does not resolve the ensemble-average image, so that |u| ^ l/0img “ we can 
approximate the visibility by expanding the exponential in 27rfu • 6: 

V{u) « {{lie)))+2niu-{{91(9)))-of 1(9))). (14) 

Here, {{f{9))) = J dfe f{9) denotes an angular average. Likewise, the visibility noise AV(u) is related to the image 
noise ^I{9) as 

AV(u)« {{AI{9)))+2Triu-{{9AI{9)))-2Tr‘^(^{u-9)'^ A/(e)Jj). (15) 

Flux Image Wander Image Distortion 

Modulation 

Thus, on short baselines, image wander will affect the visibility phase while image distortion will affect the visibility 
amplitude. Noise in |V(u)| will be quadratic in baseline: 


(lAV(u)l^) 


{{Ai{e)))^ 


1 -I- 47r^u^ 


■ {{eAije))) - 

. ■ ((A/( 0 ))) . 


(((u.0)^ AI{9)))\ 

Wm J 


(16) 


Note that the quadratic coefficient comes from the combined effects of image wander and image distortion, which act 
in opposition. The two effects could be decoupled by evaluating the second moment of the snapshot visibility, (iZg), in 
addition to the second moment of the snapshot visibility modulus (see Appendix |B.1[ ), for instance. Figurej^illustrates 
these changing contri but ions to noise in the snapshot visibility as a function of baseline length. 

By expanding Eq. in baseline, we can now easily estimate the relative strength of the position wander and 
image distortion. Again adopting a circular Gaussian source with characteristic scale rimg in the visibility domain 

(0img ~ 0.37A/riing), we And that the quadratic coefficient has a relative amplitude of « this coefficient 

is actually quite a steep function of a, ranging from 2.3 at a = 5/3 to 16.1 at a = 1.9. The aggregate rms position 
wander and distortion (in radians), y/ (A0^), along one axis is then roughly 


\/W 


A / _ 2 

nmg y 47r2(2 - a) 


{AP) (X 



^scatt 

^img 


l-al2 



(17) 
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Figure 5. Fractional rms refractive noise (Eg. |32[ | as a function of baseline for four values of linearly incremented source size, beginning 
with a point source. Left panel shows a linear scale, right shows logarithmic. The scattering parameters corresponding to these curves are 
rp = lOOro, M = 0, and a = 5/3; the sources are circular Gaussians with unscattered sizes that are approximately 0, 0.5, 1, and 1.5 times 
the scattered size of a point source. On each curve, the heavy tick mark indicates the baseline on which the ensemble average visibility falls 
to l/\/e 5^ 0.61; for a point source, this baseline corresponds to (l + M)ro. In each case, this location is near the peak of the refractive noise 
curve. The three noise regimes discussed in ^are indicated: 1) zero-baseline noise reflects refractive flux modulation; 2) short baselines 
see additional noise, increasing with the square of baseline until turning over at the noise peak, from the combined effects of image wander 
and distortion; 3) long baselines resolve the ensemble-average image and are subject to noise from scattering-induced substructure within 

the average image. An increasing source size affects the noise differently in these regimes: flux modulation falls as fractional 

image wander and distortion also as ^^/g and substructure as ^[^g- 


Note that, f or a point source, this scaling matches the result of ( [Cordes et al.|19^ (our a is their a — 2) and Romani 


et al. (1986) (our a is their /5—2). Interestingly, while most refractive effects tend to weaken with increasing wavelength, 
absolute image wander and distortion become stronger (see Table [^. 

It is perhaps more natural to express ^/{A^ as a fraction of the angular size of the ensemble-average image, 

/ \ 2 —a / \ 2—q;/2 

Then, the rms fractional wander is oc f f “ the same scaling as the flux 

modulation. But most of the decrease in fractional wander for an extended source comes from the increasing size of 
the ensemble-average image. For instance, for a Kolmogorov spectrum, the rms wander (in radians) is only suppressed 
.1/6 , 

relative to the wander of a point source. Hence, refractive image wander may remain significant even 
when the refractive flux modulation is heavily quenched. 


y’R ^img 
^scatt 


A 


mg 


by (^ 


5.3. Refractive Substructure in the Scattered Image 

Finally, observations that can resolve the scattered image will see substructure introduced by scattering. This 
substructu re co ntributes noise to long-baseline visibilities that is persistent and wideband, in accord with the conditions 
derived in (3.2 On a fixed baseline that is long enough to resolve the ensemble-average image (i.e., |b| nnig), Eq. 
shows that 


32 


(|AK(b)P)^^cx A 


oc Dd 


1 + M 

b 

1 + M 


d^y2 (y2 ■ b) |T4a ((1 + M)y 2 )\^ 


d^e |(b.V)Jea(0)r 


(18) 


This form highlights the major features of substructure on long baselines: the rms noise falls with baseline as |b|~“'^^, 
the fractional suppression of the noise due to an extended source is independent of baseline length but may depend 
on the baseline orientation, and the noise is proportional to the root mean squared gradient of the unscattered source 
brightness. The refractive noise is then inversely proportional to the squared size of the ensemble-average image. 
For the specific case of source and scattering that are circularly-symmetric Gaussians, the rms refractive noise for 

long-baseline visibilities is J (|At4s(b)|^^ °^ (^) (^) (^A^ 


From Eq. |32[ we can also quickly understand how the inner and outer scales will affect the long-baseline noise. Eor 
instance, a finite outer scale will have little effect on the refractive noise on baselines |b| ^ (1 -I- M)rA(7ut = = 

(1 -I- M)ro(rR/rout)- These baselines are long enough to resolve the decoherence length, rout, set by the outer scale. 
Likewise, a finite inner scale introduces a sharp cutoff in the noise for baselines longer than (1 -|- M)^rR. This cutoff 
arises because the scattered image consists of many magnified and distorted “subimages” of the unscattered source, 
with a characteristic scale rinj^ Thus, the visibility noise falls sharply on baselines long enough to resolve the individual 
subimages. 


® We thank Ramesh Narayan for identifying this correspondence. 
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Table 1 

Manifestations of Refractive Noise. 


Metric 


Wavelength Dependence (Point Source) 

Approximate RMS Fixed Baseline Fixed Resolution 


Fractional Flux Modulation 
Image Wander Sz Distortion 
Refractive Substructure 


2-a/2 


(roV " ( gscatt V 
V ; V ^img J 


ro 


rp 
-a/2 


il + M)roJ 


V 0img J 


2 




A^+f-s 


A final important consideration is ani sot ropic sca tter ing. Consider an anisotropic Gaussian scattering disk with an 
axial ratio A > 1. By substituting Eq. 36 into Eq. 32 we find that along the major axis of the scattering disk, the 


refractive noise is enhanced by a factor of A relative to that for a circular scattering disk with diameter equal to that 
of the major axis; likewise, along the minor a xis, the noise is reduced by a factor 1/A relative to circular scattering 
with diameter equal to that of the minor axis. ( Johnson||2013 ). The asymptotic ratio of the rms refractive noise along 
the major axis to the rms noise along the minor axis is then A“/^. 

5.4. Summary of Refractive Noise 

To summarize, visibilities on different baselines are sensitive to various manifestations of refractive noise. A zero- 
baseline interferometer sees fluctuations that correspond to flux modulation of the image. Short baselines, which do 
not resolve the scattered image, see an additional contribution: the aggregate of image wander and large-scale image 
distortion. The rms noise peaks near baselines of length rimg where 14a (nmg) = I/a/s and then transitions to a 


power-law regime on longer baselines, with the visibility noise rms falling with baseline as |AE| oc |b 


-a/2 


reflecting 


compact substructure in the scattered image. The length rimg also determines the correlation length of the visibility 
noise in the observing plane: vector baselines that differ by more than a distance of rimg will see different realizations 
of the noise. See Table for the behavior of the noise in each of these regimes. 

Because these manifestations of refractive noise are affected differently by source structure, they provide a way of 
decoupling the sour ce structure from the scattering without requiring multi-frequency measurements. For instance, 
Gwinn et al. (2014) used the level of substructure for Sgr A* at A = 1.3 cm on long baselines to constrain the pair 
of parameters I a, f/src} and found good agreement with values that have been inferred by extrapolating scaling laws 
from much longer wavelengths. Moreover, some refractive effects may be much easier to detect than others. For 
instance, refractive image wander is minimally affected by an extended source whereas flux modulation falls quite 
steeply with increasing source size. Furthermore, the effects of an inner and outer scale depend strongly on baseline. 
Thus, a comprehensive study of refractive effects can provide a robust understanding of the source and scattering 
characteristics. 


6. OBSERVATIONAL IMPLICATIONS 

Our primary result is that small-scale refractive substructure persists in the presence of an extended source, even 
on baselines that would completely resolve the unscattered source. Consequently, very-long baseline interferometry 
(VLBI) might detect structures smaller than any for the unscattered source, even if scattering does not signihcantly 
blur the source. We now give example calculations and specihc implications for two ongoing projects: RadioAstron 
and the Event Horizon Telescope. 


6.1. Implications for RadioAstron 


RadioAstron is a 10-m radio dish in a highly eccentric elliptical orbit around the Earth (Kardashev et al. 2013). 
Launched in July 2011 and operating at wavelengths ranging from 1.2 cm to 92 cm with Earth-space baselines as long 
as 370,000 km, it provides unprecedented angular resolution at these frequencies. A key science goal of RadioAstron is 
to investigate the brightness temperatures of the most compact active galactic nuclei (AGN) via direct imaging. Since 
brightness temperature Tb oc (A/0src) , the limiting brightness t emperature that a baseline can detect is proportional to 
its length squared but not to its frequency (Kovalev et al. 2005). However, the long-baseline coverage of RadioAstron is 
quite sparse because the single orbiting antenna effectively provides only a single Earth-space baseline. Measurements 
on this baseline are by far the most important contributor when inferring compact structures, so it is essential to 
quantify potential refractive substructure. We will now show that refractive noise can be mistakenly identified as a 
signature of compact intrinsic structure, leading to significant over-estimates of intrinsic brightness temperature. In 
contrast, scatter-broadening will always lead to an overestimation of the intrinsic size and, thus, an underestimation 
of the intrinsic brightness temperature. 

To derive a suitable point-of-reference, we consider a pure Kolmogorov structure function. We suppose that scattering 
material lies within the Milky Way, at a typical distance D ks 1 kpc, and that the source is at cosmological distance so 
that M = 0. A typical long RadioAstron baseline has length b ^ 10® km. We take a typical source angular size to be 
~ 300 fias at A = 6 cm, 0src ~ 1000 pas at A = 18 cm, and 0src ~ 5000 pas at A = 92 cm. These values are chosen 
to be somewhat larger than the most compact sources that RadioAstron can resolve at each frequency; the resultant 
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Figure 6. Implications of refractive substructure for the RadioAstron mission; angular broadening (top) and refractive noise {middle, 
bottom) on a 100,000-km baseline; all images utilize a Mollweide projectio n, (top) Prediction of angular broadening at A = 1 GHz for 
extragalactic objects according to the NE2001 model (Cordes & Lazio|2002l; the median angular broadening is mas x Aq^^. (middle) 
Refractive noise rms, relative to the zero-baseline flux, at A = 18 cm on a iu0,000-km baseline for a circular Gaussian source with FWHM 
of 9src = 1 mas. For comparison, the median angular broadening is 0scatt ~ 300 p.as. (bottom) Refractive noise rms, relative to the 
zero-baseline flux, at A = 6 cm on a 100,000-km baseline for a circular Gaussian source with FWHM of 9stc = 300 pas. For comparison, 
the median angular broadening is 0scatt ~ 30 pas. In each case, the rms noise varies as 0^^ ~ (^src + ^scatt) ^ • 
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9src oc A scaling is al so prototypical for self-absor bed jets ([ Blandford fc Koni^ l979) and our chosen sizes are typical of 
compac t AGN (e.g., Pushkare v fc Kovalev|2015 ). Figure [6| shows the estimated angular broadening using the NE2001 
model (jCordes & Lazio||2U02(); we take as typical values ^catt ~ 30 /xas at A = 6 cm, 0scatt ~ 300 /ras at A = 18 cm, 
and ^scatt ~ 8000 /xas at A = 92 cm, which are approximately the median of the NE2001 angular broadening and, 
thus, represent typical lines of sight away from the Galactic plane. The size of the observed, scatter-broadened image 


is 




02 

scatt' 


Under these assumptions, the fractional rms visibility from refractive noise, (Tref(b) = 
dioAstron observing bands is approximately}^ 


y^|AV;s(b)|^^, in the Ra- 


crref(b) = 0.0038 X 


= 0.0071 X 


= 0.0055 X 



I \ 5/6 

scatt \ 


-5/6 / 

105 km ) \30 ^as ) 

-5/6 / . 

/ 




img 


scatt \ 


\ 5/6 


300 fi&s 

f 06 


-2 


D 


\300 /ras ) \1000 /ras ) \1 kpc 


1 kpc 
-2 


- 1/6 


(19) 


D 


- 1/6 


-5/6 


( 0s, 


5/6 


(. 


0i, 


\8000 /tas J \ 10000 /las 


D 


1 kpc 


- 1/6 


( 20 ) 


where we have applied the relationship between rg and ffgcatt for M = 0: 

0scatt«^^^-«O.37-. 

TT rg rg 

A typical bright, compact source observed by RadioAstron might have a flux density in its compact component of 
1 — 10 Jy, so we might expect substructure with strength 5 — 50 mJy. Because baselines from RadioAstron to sensitive 
Earth antennas such as Arecibo have 4-cr detection sensitivities < 10 mJy, substructure is quite plausibly detectable, 
even on baselines that would completely res olve the unscattered s ource}^ Figure shows crref(b) as a function of 
position on the sky, using the NE2001 model (Cordes & Lazio||2002 ). 

It is evident that long baseline detections by RadioAstron must be interpreted with caution and compared with 
expected levels of refractive noise. Refractive noise would be constant over the refractive timescale of hours to weeks 
depending on the observing frequency and line of sight, and would be constant across the bandwidth of RadioAstron. 
It could be correlated among vector baselines that differ by up to a few times 10,000 km. Over many observing 
epochs, the visibility amplitudes on a fixed baseline would be drawn from a Rayleigh distribution; they would have 
a mean value of ^CTref(b) with 95% of samples drawn from the range [0.16,1.92] x CTi.ef(b). To determine whether 
a measurement of 0src is a secure indication of compact intrinsic structure, one must use Eq. [T9| together with best 
estimates of 0scatt along the particular line of sight, to see whether the measured long-baseline visibilities can be 
reproduced from refractive noise with a larger Osrc than is inferred in the absence of refractive noise. Long-baseline 
measurements that are consistently lower than the expected refractive noise could be used to infer extended structure 
that quenches refractive noise. 

6.2. Implications for the Event Horizon Telescope 

The Event H orizon Telescope (EHT ) is an ongoing international effort to establish a global 1.3-mm and 0.87-mm 
VLBl network (|Doeleman et al.||2009|). A key motivation for this development is to resolve the intrinsic structure of 
the Galactic ceiite r supermassive black hole, Sgr A*, which is significantly blurred by scattering at longer wavelengths. 
Fish et al. (2014) has shown that remaining blurring in the ensemble-average regime can be effectively removed at 
A = 1.3 rnrm However, while angular broadening decreases steeply with increasing frequency, the fractional effect 
of refractive noise becomes stronger with increasing frequency, so the influence of refractive substructure remains an 
important consideration for the EHT. Moreover, because Sgr A* only transitions to weak scattering (rg = rp) at 
^2 THz, refractive substruct ure will affect all foreseeable VLBl of Sgr A*. 

Indeed, Gwinn et al. (20141 detected refractive substructure in the scattered image of Sgr A* at A = 1.3 cm, apparent 
as persistent visibility 3n baselines that resolved the ensemble-average scattered image. The most pronounced effects of 
substructure will occur when the angular broadening is comparable to the intrinsic structure (for Sgr A*, at A ^ 5 mm). 
Figure shows the expected refractive noise for Sgr A* at wavelengths and baselines appropriate for the EHT. The 
noise is si g:nificantly affected by the intrinsic sour ce structure, which is currently only well-constrained in the East-West 
direction ([Doeleman et al.|2008| |Fish et al.|2011|). With the assumption of a 40 /xas « 4i?gch circular Gaussian source, 
the refractive noise at A = 1.3 mm is ~60 mJy for long East-West baselines and is ~25 mJy for long North-South 
baselines (the difference is because of ~2:1 anisotropy of the scattering). Continued observations are essential to better 
characterize these noise properties. 

Note that while current EHT measurements on long baselines are similar to the level of refractive noise expected 
for a point source, they cannot be entirely the result of refractive substructure because the ensemble-average visibility 

® At A = 1.2 cm, the median angular broadening predicted by the NE2001 model is only pas, resulting in rg ~ rp for a scre en a t 
D = 1 kpc. Thus, most observations with RadioAstron at A = 1.2 cm will be in a strong/weak scattering transition regime and Eq. |19| is 
not directly applicable. 

^ For estimates of the RadioAst ron SEED in each obs erving band, see the RadioAstron User Handbook (www.asc.rssi.ru/radioastron/ 
documents/rauh/en/rauh.pdf) or [Kovalev et ah] ([2014^. 
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Figure 7. Noise from refractive substructure expected for Sgr A* at A = 1.3 mm {left) and at A = 0.87 mm {right). For each wavelength, 
we show the expected noise for both a point source and for a 40 /ras (FWHM) circular Gaussian source. Estimates given in mjy assume 
a 2.4 Jy source in both cases. At the shorter wavelength, A = 0.87 mm, the refractive noise for a point source is stronger but is also more 
heavily quenched by the 40 /ras source. Note that these estimates do not change sharply with different scattering assumptions, such as 
a break to Kolmogorov scaling at A~1 cm because the increase in refractive noise from somewhat weaker scattering is partially offset by 
the concurrent increased quenching from source structure. However, these estimates demonstrate that the refractive noise does depend 
sensitively on the intrinsic source structure at these wavelengths. 


would then be far higher than the measured values (i.e., the intrinsic structure would need to be more compact than 
has been inferred). The consistency of the EHT measurements over observations in different years further argues 
against refractive noise as a dominant contributor. Nevertheless, because the estimates shown in Figure anticipate 
refractive noise at A = 1.3 mm that is higher than expected thermal noise in EHT measurements (see, e.g., |Lu et al. 


2014|), refractive noise will be an important consideration for the interpretation of forthcoming EHT data, especially 
on baselines for which the ensemble-average visibility is low. Although this noise would remain constant over a single 
EHT observation and across the bandwidth of the EHT, by co mbining different observations one could reduce the 


20141. 


noise and approach the ensemble-average limit (Fish et al. 

Refractive imag e distortion may al so cause distortion of the nearly circular photon ring that surrounds the “shadow” 
of the black hole (Falcke et al.||2000[), which rnay complicate tests of general relativity based on potential asymmetry 


of the photon ring (e.g., Johannsen &: Psaltis 2010|). While refractive distortion would be stochastic and would not 
change the shape of the shadow when averaged over many observing nights, it may be an important consideration for 
single-epoch studies. The supermassive black hole in M87, which subtends a similar angular size but has negligible 
effects from scattering at EHT wavelengths, will provide an important point of comparison. 

Improved imaging strategies may help to mitigate the effects of refractive noise. However, recall that the correlation 
length of the noise intro duced by substructure is approximately the location at which the ensemble-average visibility 
falls to l/^/e (see §3^1. For a 40 ftas circular Gaussian intrinsic source, this length is approximately 2200 km at 
A = 1.3 mm and 1600 km at A = 0.87 mm. Because the completed EHT will have E ast-West baseline s extending to 
^8000 km and North-South baselines extending to ~11,000 km (see, e.g.. Figure 1 in Fish et ar]|2014 1, there may be 
many uncorrelated elements of the noise in a single observation, complicating mitigation strategies. 


7. SUMMARY 

Refractive scattering poses a complex observational challenge. It introduces effects that are long-lived and wide¬ 
band, and that can be difficult to unambiguously identify without multi-epoch observations. We have explored the 
contributions of a particul ar refrac tive effect: scattering-induced substructure in images. Building upon the theory and 
methodology developed by |GN89l we have shown that small-scale substructure persists in the presence of an extended 
source. In particular, the action of scattering is not a convolution in the average-image regime, which describes nearly 
all foreseeable VLBI of extended sources. As a result, when the scatter-broadening is comparable to or exceeds the 
intrinsic structure, the signature consequence of refractive substructure - elevated visibility on long baselines - can 
persist even on interferometric baselines that are long enough to resolve the unscattered source. These long-baseline 
visibilities are not an indication of compact source structure. This result starkly contrasts with the behavior of the 
diffraction pattern in the observing plane, which is smoothed by an extended source by convolution. 

Our results have immediate implications for VLBI studies of galactic nuclei, including studies of AGN with RadioAs- 
tron and of Sgr A* with the EHT. Refractive noise is an especially important and subtle consideration for resolved 
images when the scatter-broadening is comparable or somewhat subdominant to the intrinsic source structure. The 
refractive noise limits direct imaging to a maximal resolution determined by the scattering and can introduce spurious 
compact features into images. Nevertheless, our results also indicate that these instruments may provide valuable 
insight into the turbulent ISM from the new perspective that refractive substructure affords. For instance, substruc¬ 
ture reflects properties of the turbulence on much larger scales than those that produce the angular broadening. Also, 
unlike refractive flux modulation, refractive noise on long baselines does not need to be disentangled from intrinsic 
source variability; it may, therefore, be a more robust observable than flux modulation and an easier measurement 
than position wander (which requires absolute phase calibration) or image distortion (which requires precise ampli- 
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tude calibration). Refractive noise can also decouple the intrinsic s ource structure from scatter broadening without 
requiring assumptions about the frequency dependence of each (e.g., Gwinn et al. [2014 1. 

Ultimately, even refractive noise may not determine a fundamental resolution limit tor VLBI. Famously, scattering 
processes do not destroy information, but rather they add information until the original information is obscured. 
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APPENDIX 

A. THE PHASE STRUCTURE EUNCTION 


To encompass the phase structure function in a simple functional form, we use the following expression (GN89 
Eq. 3.1.1): 


D^iy) 





if |r| < Tout 
if |r| » rout- 


( 21 ) 


In this expression, C (X v , reflecting the behavior of the index of refraction for frequencies much greater than 
the plasma frequency (kHz in the interstellar medium). Thus, for a pure power-law, the phase coherence length 
{D^{ro) = 1) satisfies tq oc Eor a negligible inner scale, rin ^ rg, we have 


On the other hand, if r-m ^ tq then 




(S)’ 

(t) (S)‘ 


if r < Tin, 

if r » Tin. 


( 22 ) 


(23) 


While the power-law form at large r is equivalent in the two cases, the correspondence between the C and rg is affected 
by an inner scale. This change is what causes the minor modifications to refractive noise on baselines shorter than 
^ from an inner scale. 

Einally, we account for anisotropic scattering for a pure power law using the form 

/ 2 2 

i3*(K!/})= V + Y ■ (24) 

\'0,x '0,y) 

B. DETAILED GALCULATION OF REFRACTIVE NOISE 

We now calculate the variance of the refractive noise (i.e., the noise At4 = 14 — (14)ea = 14 — 14a in the average 
visibilit y). Th is variance is used in to quantify various manifestations of refractive noise. Our treatment closely 
follows GN89| although we do not restrict ourselves to long baselines or isotropic scattering, and we include effects 
from an exten ded s ource. _ 

We first, in §B.1[ derive a general expression for the second moment of the snapshot visibility. Then, in §B.2[ we 
isolate the contribution of refraction noise and obtain a simplified expression for it in the strong-scattering regime. 

B.l. The Second Moment of |Uss| 

To calculate the variance of the snapshot visibility modulus, we must multiply the expression for the snapshot 
visibility (Eq. by its conjugate before evaluating the ensemble average over the screen phases. The ensemble- 
average over screen phases only affects one term in the resulting integral, which we can readily reduce, as for the 
ensemble-average visibility (j |3.3| ): 

gi[<^(xi)-<^(x2)-<^(x3)-|-0(x4)] \ _ g- 5((0(X1)-0(X2)-<^(X3)-|-0(X4))^)^^ 

/ ea 

_ 5 [-D0(Axi3) + D^(Ax24)+-D0(Axi 2) —i^0(Axi4) —D0 (Ax23) + D^(Ax34)] 


(25) 
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Figure 8. Refractive and diffractive contributions to the second moment of the snapshot visibility. 


where — Xj. The second equality follows from 

(</)! - ^2 - ^3 + ^ 4 )^ = (</'! - + (</>2 - + ih - 4'2f - {(t>l - (t>if - (^2 “ </> 3 )^ + (<>^’3 “ ( 26 ) 

Using a Hadamard transform, we then change variables to exploit this pairwise representation: 

a 1 1 1 

1 i' -1 :1 ■ 

1 -1 -1 1 

Finally, we can apply the identity x\ — x\ — x\ + x\ = 2 (y;^ • y 4 + y 2 ' Ys) to obtain 
\ / ea (27rrp) J 

^ g -5 [-D0(Axi3)+D,^(Ax24) + -D(6(Axi2)--D(6(Axi4)-£I(6(Ax23)+-D4.(Ax34)] 

X Usrc(- (1 + M) (y2 + y 4 ))v;;^((l + M) (y4 - y2)). 


/yA 

/ 

y 2 


y3 

\y4/ 




(28) 


Because the Ax^- can be written using only {y 2 3 4 }, we can integrate Eq. 28 over y^ to replace the exponential in y^ 
by (2'Kr‘^Y We then integrate over y 4 to leave 


(|v;,(b)p) = - 2 / ^^2 3 'b]-y3e-5[2i54y3)+2i?Ry2)-0Ry2+y3)-0Ry2-y3)] ((1 + M)y^)\^ . 

( 27 rr|) J 

(29) 


Note that the phase structure function is still completely general, as is the intrinsic source structure. Moreover, we 
have not yet made any approximations based on the strength of the phase fluctuations. 


B.2. Refractive Noise 

Three effects contribute to a snapshot visibility: diffractive noise, refractive noise, and Vea{b). However, for the 
second moment, as expressed in Eq. |29[ each arises from a different region of the integral. In the strong scattering 
regime with a < 2 , the dominant contributions to the integral have at least one of {|y 2 l i lysl} ^ I’o- The region 
{|y 2 l ^ 0 : lysl ^ 'I’o} gives the diffractive noise, as is evident by its suppression from an extended source. The 

remaining region of the integral gives the second moment of |I4|- Figure ^illustrates this geometry. 

To isolate the noise At4) note that the contribution of I4a can be written as 


(27rr2)' 


J d^y,^,N^rny. + (l+M)-^b].y,^-D,iy,) {{1 + M)y^)f = |14rc (b)|' = |Vea(b)|". 


Because they focused on long baseline properties, (GN89I neglected the phase variations 6 *’'^ X 2 y 3 in Eq. 


29 


(30) 


which 


are then small over the region that dominates the refractive contribution, giving an 0((1+ M)rn/b)^ corr ection. With 
this approximation, and in the strong scattering regime, the integrals in Eq. 30 factorize ( Johnson||2013 ), simplifying 
the remaining calculation. However, we must retain the coupling term in order to understand the retractive noise on 
all baselines. 
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Johnson & Gwinn 


The integral over yg is the one that is challenging numerically, but we can get it into a manageable form. Because 
D^{y) = Dcj,{—y) and [yg] |y2| for the average visibility, we have that 


2 £'<>(y 3 ) - D^{y2 + ys) - ^0(y2 - ya) ~ - (y2 • D<i>{y)\ 


y=y3 


(31) 


Importantly, this rep rese ntat ion does not assume that D^{y) is isotropic. 
Combining Eq. Eq. and Eq. then gives 

(|AK(b)P)^^=(|l(,(b)n^^-|Eea|^ 


( 27 rr|) 


1 


( 27 rr|) J 

(27rr|)^ J 
(27rr|)^ J 

where we have introduced 


d y2,3e''F 




A [5^2+ i+M]'y3 


^-l[2D^(y^)-D^{y^+ys)-D^{y^-y^)] _ G g-D^y^) ((l M)y2)| 


[y2+TTTu-]'y3 


i (y2 • V)^Zl0(y) 


y=y3. 


|Eea((l + M)y2)|' 


C«V2,3 ^ (y 2 • [y 2 + (1 + ^) ^"Zl,^(y 3 ) |v;a((l+M)y 2 )|^ 

d‘^y2 ^ (y2 • [y2 + (1 + M)-^h]Y {y^ + (1 + M)-^h) |Eea ((1 + M)y2)|^ , (32) 


D'j>iy)= I c^V3e^^^"£’0(y3)- 


(33) 


Once Drj,{y) is specified (analytically), the integral in Eq. 32 can easily be performed numerically. 

Note that DAy-i) is closely related to the power spectrum of phase fluctuations a t a scale q = yg/r| (see, e.g., 
Tatarskii||1971[ Goodman fc Narayan||19^ Coles et al.||1987 Lambert fc Rickett]|2000 1. Specihcally, after removing a 
constant phase offset, which does not affect the snapshot visibilities, (yg) oc —Q (q = yg/r|) . Hence, positivity of 
the power spectrum enforces positivity of the refractive noise. For instance, consider an isotropic power-law structure 

function: D^{y^) = ^ . Then 


n z- nl-l-Q r(l-|-Q:/2) 

^ (^3) - -2 ™r(l-a/2)''o 


^3 

„2 

'F 


-{2+a) 


We can easily generalize this result to include inner and outer scales in or anisotropy. For instance. 


^ T{-a/2) ° 


-k 




_(l+a/2) 


jn ZS 

^0 I 


\ ^out/ 

Similarly, to account for anisotropic scattering, we use the form for D^{{x,y}) given in Eq. 


24 


Then 


D^{y) = -22+“7r 


r(l-ka/2) 


{rl^x^ + rlyv"^) 


2. -(l+a/2) 


r(-a/2) 

C. SCATTERING SIMULATIONS 


(34) 


(35) 


(36) 


have been performed in many contexts (e.g., Coles et al. 

1995 


2010 

Habibi et al. 

most similar to the one-dimensional simulations of|i\G89[ with tJ 

le primary dilterence 


Figu re [3| Similar simulations 
2^, but ours are perhaps 
:hat we are simulating a two- 

dimensional scattering screen. 

We first generate an Ax A grid of independent zero-mean, complex Gaussian random variables. Next, we imprint the 
appropriate power spectrum Q of phase fluctuations by taking the discrete Fourier transform of the phases, multiplying 
by VQ(q)) and then inverse Fourier transforming the result. Finally, we calculate the empirical structure function 
and apply an overall normalization to the phases so that D^ir^) = 1. 

For the power spectrum, we employ the commonly used form 


<3(q) fX (|qp +9min) (omix) , 


(37) 


which corresponds to a phase structure function with an index a, inner scale rin ~ l/(Zmax, and outer scale rout ~ l/^min- 
Note that specific conventions for the form of Eq. 37 and the relationship between rju/out and (7max/min vary throughout 
the literature. 
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After calculating the screen phases, we can evaluate the electric field in the observing plane (Eq. f^, snapshot 
visibilities (Eq. [^, or the scattered image (Eq. [Io| ). For the scattered images shown in Figure]^ we cmculated the 
phase gradient at eac h lo cation using a discrete approximation with nearest-neighbor grid points and then utilized the 
approximation of Eq. |12| 

The greatest computational hurdle is then storing the random phase screen in memory; the number of screen phases 
that must be stored is ~ (i"R/ro)^ = (rp/ro)^. To achieve the largest ratios rp/rp 3> 1, appropriate for strong 
scattering, we set rg to be equal to the grid spacing. Our simulations have a linear size of 2^^ points, so we can readily 
achieve rp/rg 35. However, because the phase varies linearly on scales shorter than ri„, one only needs to resolve 
the larger of {rg, ri„}. Hence, simulations could be performed with a much larger rp/rg if they also included a large 
inner scale. 
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